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Abstract 

Explicit general formulae for the tensor reduction of two-loop massive vacuum 
diagrams are presented. The problem of calculating the corresponding coefficients 
is shown to be equivalent to the problem of constructing differential operators (pro- 
jectors) extracting the coefficients of the momentum expansion of massive scalar 
three-point functions (with any number of loops), so the solution to the latter prob- 
lem is also given. 
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1. Introduction 

Tensor reduction of Feynman integrals containing loop momenta with uncontracted 
Lorentz indices in the numerator is very important for various realistic calculations in the 
Standard Model (and beyond). For one-loop diagrams with different numbers of external 
lines, several approaches and algorithms were developed [|I], 0]. For two-loop self-energy 
diagrams, the problem was considered e.g. in ref. 0, whilst the three-point two-loop case 
is more complicated and requires results for the integrals with irreducible numerators [Q. 

The problem of finding general algorithms for the calculation of two-loop vacuum 
diagrams is interesting because these diagrams occur in many important physical prob- 
lems. We can mention, for example, calculation of the two-loop effective potential in the 
Standard Model |^ and two- loop contributions to the p-parameter P, 0. Furthermore, 
consideration of many other problems requires constructing efficient algorithms for the 
momentum expansion of two- loop massive diagrams with two 0, three 0, or more ex- 
ternal legs. The coefficients of such expansions involve massive vacuum diagrams with 
tensor numerators. A procedure for reducing the tensors to scalar numerators enables 
us to calculate these coefficients analytically, since general results (including the relevant 
terms of the e-expansion) are known for the corresponding scalar diagrams P, ^, ^, |10[ . In 
the two-point case, the reduction is relatively simple because, as all indices are ultimately 
contracted with the same external momentum, only completely symmetric tensors need 
to be considered. However, starting from the three point-case, when there are two or 
more independent external momenta to contract with, a complete tensor decomposition 
is required. 

The answer is usually easy to obtain when one has only a small number of integration 
momenta in the numerator. The situation becomes more complicated when one needs 
to calculate the higher-order coefficients of the expansion. This is the case when we are 
approaching the physical singularity (threshold), where the expansion still works but the 
parameter we expand in is already not so small. Moreover, in the approach (which 
involves conformal mapping and Fade approximations) , in order to get accurate numerical 
results beyond the threshold(s), one needs to calculate expansion coefficients of very high 
order, and therefore tensor integrals of very high rank. These are some of the reasons 
why general algorithms for the tensor decomposition of two-loop vacuum diagrams are 
important. 



We would like to mention recent progress in this direction. In ref. |Tl|, a closed ex- 
pression, in terms of gamma functions, was given for arbitrary tensor vacuum integrals in 
which two of the masses are equal and the third is zero. Furthermore, it was noted that 
the algorithm used in that paper for reducing the tensors to scalars can be extended to 
more general cases. Some relevant formulae for diagrams with arbitrary masses, where the 



numerator is contracted with one or two external vectors, can be found in refs. 1121 , [131 , [14 
Uncontracted tensor integrals were considered in refs. [|1^, [1^ . In , a general expression 
for the matrix inverse to one containing the coefficients required for the tensor decompo- 
sition was constructed, while in [[T3| a version of the direct formula was presented. While 
solving the problem in principle, these results |1^, [1^ could not completely satisfy us 
because they were not explicitly symmetric and the coefficients required for the decompo- 
sition were rather cumbersome, so that their use in computer programs was less efficient 
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than recursive calculation of the corresponding coefficients. In this paper, we present a 
simpler general solution to this problem. 

Our results also provide a complete and explicit answer to a different question, which 
is related to an alternative approach to the expansion of scalar three-point diagrams. 
Instead of Taylor-expanding all propagators, one applies certain differential operators to 
the diagram and then sets the external momenta to zero, in order to obtain the Taylor 
coefficients of the diagram as a whole. In this way, one does not need to evaluate any 
tensor integrals, but is faced with the problem of finding the correct differential operators 
to use. We shall show that this problem is, in fact, equivalent to the tensor decomposition 
of two-loop vacuum integrals, which enables us to give a closed expression for the operator 
that projects out the Taylor coefficient of any given order in the three external kinematical 
invariants. 

The remainder of this paper is organized as follows. In section 2 the notation is 
introduced and the connection between two-loop tensor reduction and projectors for three- 
point functions is established. The solutions to both problems involve a set of universal 
coefficients, which are calculated in section 3. A summary and a discussion of the results 
are contained in section 4. 

2. Definitions 

In this paper, we shall use the following notation: 

/ [something] = J J d"p d"g {something} F (j)^ , q'^ , (pq)^ , (1) 
and we are interested in expressing the tensor integrals 

PfJ.1 ■ ■ ■ PlJ.Ni ^0"l ■ • ■ 9'o-jV2 

in terms of scalar integrals. In eq. (^, n is the space-time dimension |jl6| and F{p'^, g^, (pq)) 
is an arbitrary scalar function depending on Lorentz invariants of the loop momenta p 
and q. Usually, it is a product of (powers of) propagators, 

(p'^-ml) [q^-ml^ ''^ (^{p-qY-mfj , (3) 

times a polynomial in p^, and {pq), but the concrete form of this function will not be 
important for our discussion. 

Because F is an even function of p and g, it is clear that the integral (|^) vanishes when 
Ni + N2 is odd. When Ni + N2 is even (from now on, we shall always assume this), since 
we have no external momenta, the result must be a Lorentz invariant tensor made out 
of metric tensors, and it must be symmetric in the two subsets of indices, {/ii, . . . ,^Ni} 
and {(Ji, . . . , (Tatj}. A basis of tensor structures with these properties can be described as 
follows. 

Each structure is characterized by three integers, ji, j2 and js, such that 2ji + J3 = A^i 
and2j2+j3 = N2. It can be constructed by taking a product of ji metric tensors Sf^i^^-, J2 
tensors g^.a^ and j's tensors S'/i^o-j., and then symmetrizing it in the /x's and a's by taking 
the sum of all distinct products of metric tensors we get from it through permutations of 
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the /i's and of the ex's: 

+ permutations . (4) 
The number of terms on the r.h.s. is 

^^^^^^ - 2n^n j^lj^lj^l ■ 
Equivalent representations of the tensor structures are 

= h,, . . . k,^^^ fc^., . . . k,^^^ , (6) 

where we have adopted the notation of [0| for d'Alembertian operators, 

For simphcity, we have suppressed the Lorentz indices on the Lh.s. of (|^). We note that 
in rT|^ d'Alembertians were used in a similar way to generate the required tensor 
structures. 

At given A^^i and A^2; the number of independent tensor structures (|) is 

T{Ni, N2) = l + min {[Ni/2], [N2/2]) , (8) 

where [Ni/2] is the integer part of Ni/2, and each of these structures is already determined 
by one of the j's, e.g., j^. Because we shall often need to sum over the set of all tensor 
structures, we introduce the following notation: 

E - E ^ (9) 

{j} jlj2,i3 

2ji+j3=Ni 

emphasizing that this is actually just a single sum. 

Now, we return to the integral (0) and write it as a linear combination 

-' PfJ-l ■ ■ ■ PfJ.Ni 9(71 ■ ■ ■ Q'o-JVj 

(10) 

{J} 

with some scalar coefficients Ijijijs which are to be determined. By contracting (PUj ) with 
each of the structures (^) we obtain the following system of T{Ni,N2) linear equations 
for the /jjj2j3's: 



^jlj2j3 ^ 



{f} 



Here, we would like to mention that eq. (B.4) of |0 is equivalent to the result presented in if 
one explicitly calculates the result of applying the d'Alembertians. 
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where 

It is useful to think of the coefficients Xjij2h;j[j'2j'3 elements of a T{Ni,N2) x 

T(Ni,N2) "contraction matrix", whose rows and columns are labelled by the triplets 
(ji)j2, is) and (ii, ^2)^3)) respectively. We shall denote the inverse of this matrix by 0, 
i.e. 

b"} 

Since x is a symmetric matrix, is also symmetric: 



(14) 



Using to solve the system of equations (|n|), we can re- write eq. ( jTDD as 



{j} m 

(15) 

Thus, we have reduced the tensor integrals (^ to a combination of scalar integrals (car- 
rying the same total powers of momenta p and q) multiplied by the tensor structures 
{ji, j2, js}- In the next section, we shall derive explicit expressions for both and x- 

Here, we address the second question posed in the introduction. Consider a massive, 
scalar, three-point Feynman diagram depending on two independent external momenta 
ki and k2, and suppose we are interested in the coefficients of its Taylor expansion, 

c{k,,k2)= E c,,,2rAkiy'ikiy%hk2y^ . (le) 

jl j'2 J3 

Our goal is to construct a set of projection operators that give us the coefficients Cj^^j^j^, 
using derivatives with respect to the vectors ki and k2. This will allow us to apply them 
to the integrand of the Feynman integral through which C{ki, ^2) is defined. 

Following let us apply a product of d'Alembertians to both sides of (|16D, and then 
set fci = ^2 = 0. This gives 



□fi nil C{k,, k2)l^^^^^^ = E 9;.^.^°n nil {kly'^{klyKkik2y^ 

^ {f} 



nl n' 



mm.Y. f 0,^^, , (17) 



where iVi = 2ji + and N2 = 2j2 + js- The second line follows from eqs. (|^) and (|T^). 



Solving the system of equations (|T^ gives us the operators we are looking for: 



c,u2.3 = ^ffi E 0.U2.3;.;.^.^ Hj^,:, °il °2l nil 



fcl=fc2=0 ' 



where the coefficients 0jij2j3;jij2J3 same as the ones needed in the tensor reduction 

formula (^)! For the coefficients Coojg and CjjOjs; this result, combined with eq. 
presented below, coincides with the expressions given in H. 



Another way to see the connection between (|15[) and (|I8[) is by taking the scalar 
function C{ki, /C2) in the latter to be / (kip)'^^ {k2q)'^^ , and noting that 



P^ll ■ ■ ■ Pt^N^ 1' 



1 



Ni\N2\ dk^' ...dk^' dkl 



{hpf'{k2q) 



N2 



(19) 

From this point of view, it is also clear that the connection can immediately be generalized: 
the 0's occurring in the tensor reduction of L-loop vacuum diagrams are the same as in 
the projectors for the Taylor coefficients of {L + l)-point scalar diagrams. 

3. Results 

We calculate the coefficients ^ju^js^iiaia inductive method, based on recurrence 

relations, which enables us to express them in a compact and symmetric way. Afterwards, 
we shall obtain x by inverting 0. 

We start from the following simple formulae for contracting two of the Lorentz indices 
of the tensor structures (Q), which can be derived, e.g., by using (||): 

9m^a^^{jld2,3?.] = (ra + 2ji + 2j2+j3-l) {jl,j2,j3-l} + (j3 + l){jl-l,j2-l,j3 + l}, 

(20) 

9^.^,^^^^M, {JUJ2,J3} = (n + 2ji + 2j3-2) J2, J3} + 2 (J2 + I) {ji, J2 + I, J3-2} (21) 

and an analogous formula for contraction with 5'o-^2-io"jv2 • r.h.s., the structures 

{ji)j2, is} should be set to zero whenever any of the j's becomes negative. Using these 
contractions, the following recurrence relations for ^jijaiayjjajg '^^^ be obtained: 

A^iA^2 {in + Ni + N2- - 1) (l^nMUl^j!, + Us - 1) '/'iiiaisiiJ+U^+i^-s} 

= J3 0iij2,j3-i;i; ifs > 0) , (22) 



ATi (iVi-l){(n + Ari+j^-2)0,, 



= 2 Ji 0ii-ij2,i3;j(-ij^4 (^'^ > 0) ' (23) 



and also a relation similar to ( [2B| ) but with interchanged indices, 1^2. If either the 
second (p on the l.h.s. or the one on the r.h.s. has a negative index, it is to be set to zero. 

Let us consider some fixed values for ji, j2 and and suppose we already know the 
0's on the r.h.s. of, say, (p2D, corresponding to (A^^i — 1) /i's and (A'^2 ~ 1) cr's. The number 
of equations (|2^) is T{Ni — 1,N2 — 1). If the number of tensor structures does not increase 
between (iVi — 1,A^2 — 1) and {Ni,N2) (i.e., when A^i and A''2 are odd), this means we 
have enough relations to determine the 0's on the l.h.s. completely. On the other hand, 
if at this step one extra tensor appears, we need one more relation between the 0's. The 
same is true if we want to use (^) to go from (A^i — 2, A^2) to {Ni, N2). The additional 
information can be obtained by contracting (|r^) with A;^^ . . . /c^^i k'^^ . . . k"^^ and using 
(cf. eq. (B.IO) of Q) 



{kpf^kq) 



N2 



{k 



2\{Ni+N2)/2 



2{Ni+N2)/2(n\ 

\2j(Ni+N2)/2 



3Z 



(24) 
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where {a)j = r(a + j)/r(a) is the Pochhammer symboL Comparing the coefficients of 
I[{py^{qy^{pqy^], we see that 



The recurrence relations (p^)-(]23|) can be solved in several ways, leading to different 
representations for 0. First, we consider the special case when ji = 0, so that the r.h.s. 
of eq. ( p3D vanishes and we get a relation between ^oj^iayJiaJs 'i^oj2i3;j(-i,i2-iJ3+2- 
Repeated application of this relation yields 

J2!(— +Jl+J3j., 

To determine the normalization, we insert this result in (p5D . The sum is a terminating 2-F1 
series of unit argument, which can be written in terms of Pochhammer symbols. Substi- 



tuting the expression for 0oj2j3;Oi2j3 ^^us found back into (|26D and using the symmetry ( p!4D 
gives 

'Phj2nMN2-m)/2,m - / s • (27) 

To handle the general case, we move the second term in (|23|) to the r.h.s. In this way, we 
get a recurrence relation in j[, which has the following solution: 



Z=max(^Oji-ji j 

T^X (Pji~j[+i,j2,j3;0,j'2-i,j'3+2i ■ l^oJ 



2 V2yjj+j^+i 



We can write this as a familiar sum over j's, (P), by defining = j[ — /, J2 = J2 ~ ^ ^'^'^ 
is ~ is + 2/. Substituting the result for j[ = (P7|), we find 



^ ^ Ai3+i^')/2V 2 ^29) 



In (|29|) and in the rest of this section, it is understood that terms with factorials of 



negative integers in the denominator vanish. Note that (EOf) clearly shows the symmetries 



3i32rivj[i'2i'i V i[3'^3'y3-Lj2ri r'j^jusyailis" 
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Another special case for which we find a simple expression is when = 0. This time, 
the r.h.s. of (p2|) vanishes, and steps similar to the ones described above lead us to 



(f)'(f)'(-ir^'^^3!(^)„. 



+N2-j3)/2 



(30) 



{Ni+N2)/2 



Using (p2D as a recurrence relation in jg gives the following result: 



(=maxlOji-jJ 



n-1 
2 



^ /.„_i\ ' ' . (f'ji,j2,j3-j'3+2l;j[+ld'2+l,0 ■ (31) 



2 ^"'^2j(+2i^+j-+M 

By writing j[' = j( + /, = jj. + ^ ^'^d is = Js ~ 2/, and inserting (pOD we obtain a second 
representation for 0: 



-i"/ /i'/ -i"/ 



E 



^ V 2 Jjn+j. \ 2 ; (^■^+j,+j-»+j»)/2 V 2 ;(j/+jV+j»+j»)/2 ^^^^ 

The representations (|29|) and (|32D complement each other. The number of non- 



vanishing terms in (poD is 1 + min(ji, j2, jj, ^2)5 whereas the number in (|32D is 1 + 
min([j3/2], [j3/2]). Therefore, depending on the values of the j's and j"s, it can be 
much more efficient to use one or the other, even though both lead to identical results^. 
A simple example is when A^^i and A^2 are odd and jg = 1. In this case, (|3^) consists of 
just one term, with jg = 1, which gives 

, _ { 2 )■{ 2 J3-[ 2 ;(^^+jv,_,3_i)/2 

<Pjij2j3;{Nj-l)/2,{N2-l)/2,l 



V 2 y V 2 J[Ni-l)/2\ 2 J {N2-'1)/2\^J (Ni+N2)/2 

(33) 

An interesting feature of (p9|) and (p2D is that, in each term of the sum, the dependence 
on j and on j' is factorized. Both formulae have the structure = M'^DM, where M is 
a triangular matrix and D is diagonal. The matrix elements of the matrices M, and of 
their inverses M~^, can be presented as 

_ {a)j+i , .1. _ (-1)^-' (g + l)2j 1 . . 

[l-jy. [i-jy- (a)2i (a + l)i+« bj 



4 



It is possible to verify that (29) reduces to (|3C|) when 73 = 0, and ( p2| ) to (27) when j( = 0, with the 



help of a theorem on the summation of 5F4 series; see §4.3, eq. 3, or Il9l], app. Ill, eq. 13. 
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With this decomposition, it is now easy to invert the matrix 0, which gives us the following 
representations for the contraction matrix x'- 

(35) 

if we use (Eol), and 



(36) 



if we use (^2p. We also note that the determinant of is given by the product of the 
diagonal elements of M^, D and M, 

det0 = n ^ , ^ 7% . (37) 

\hNi\ N2\ (n-2), f^) f^) 

The matrices M have a rather nice interpretation. They transform the basis of tensor 
structures {ji, j2, js} into orthogonal bases. The orthogonal tensor structures involved in 
our first representation, (^) and (|35D, can be written as 



(^)! 



(^^2^)iVl/2(^^2^|7V2/2^{n-2)/2 / (^1^2) 



where 



. . . 9A:r^ dk^,^ . . . dk^ ' ' " \x/klkl 

^^./2 ^.2/2 ^(„-2)/2 I _ _ ^ ^ (3g) 

°11°22 

[i/2] f-iVr^i 

^7(-) = E (2-)^'"" (39) 

are Gegenbauer polynomials. If we remember that A^^i = 2ji + and N2 = 2j2 + js, we 
see that (pSD does not really contain negative powers or square roots of d'Alembertians. 
The orthogonality of the structures (^) can be explained using the fact that 



n,,{kiy^/\kiy^/'c\:^''^^' f ] = o (4o) 



(and the same for □22)? and it is closely related to the orthogonality of the Gegenbauer 
polynomials with respect to integration over the unit sphere in n-dimensional Euclidean 
space. 
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An analogous explanation can be given for our second representation, (|3^) and (pq). 
The corresponding orthogonal tensor structures are connected with polynomials that are 
annihilated by □12, 



1=0 



l\{ji-l)\{j2-l)\ \ klkl 



(41) 

For completeness, we mention that a third representation for ^ju^iayjjais "^^^ ^le derived 
by using (|23|) in the opposite direction, as a recurrence relation in jg. It is a sum of 1 + 
min(ji, [j3/2]) terms, in which the dependence on j and j' is factorized. However, it is less 
symmetric than (|29| ) and (|32D and does not appear to have such a simple interpretation. 
It is also possible to write ( ^2]) as a recurrence relation in j[, but it is more difficult to 
solve than the other ones. 

To conclude this section, let us insert our first expression for ( pQ]) into the reduction 
formula (pIS]). The result can be written in the following wayQ: 



PfJ.JSl, QcTl •••<?( 



1 • ■ ■ l^f^Ni 



1 



E 



J3! 



J3 



X / 



(-2).3(f),,,Jt) 
{{iVi,iV2;j3}}- 



i2+j3 



(42) 



When contracted with ki and A;2, this gives 



N2 



Nim^ is! 



E 



(f) 



J3 



J3 



n - 2),. 



X J 



(^2)^,/2(^2)^./2^(n-2)/2 ('Jpi_')j (^2);V./2 (^2) ^./2^(n-2)/2 | 



2i,2 
2, 



and, similarly, the three-point functions ( [T^ ) can be expressed as 



C(fci,A;2) 



E 



E 



J3! 



32+33 



X 



□ 



Afi/2 |--|Af2/2 ^(n-2)/2 



^-2^2 
"'1"'2 



11 



22 



□ 



12 



J3 



°11°22, 



C{kr,k2 



ki=k2=0 



(43) 



(44) 



^ In eq.(6) of ref. |T^, the scalar coefficients are also written in terms of Gegenbauer polynomials, 
but a different basis is used for the tensor structures. By some transformations, and the summation of a 
hypergeometric series of unit argument, one can show that the result is equivalent to p^). 
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4. Conclusion 

The main results of this paper are the explicit formulae (p9D and (|32D for the coefficients 
%ii2j3;j(j2J3' "^hich necessary for both the reduction of two-loop vacuum tensor integrals 



of arbitrary rank to scalar integrals ([T5|) , and for the construction of the projectors (|n 
that extract the Taylor coefficients of scalar three-point diagrams. The tensor integrals 
are written in a more compact way (^21) by using orthogonal bases involving Gegenbauer 
polynomials for the tensor structures and scalar integrals. A similar representation is also 
derived for the three-point functions, eq. (|4^). This shows, once again, that when one 



expands in {kik2), the basis of the Gegenbauer polynomials is useful PU| (see also in ||TB|) 



On the way, we have also obtained explicit expressions for the elements of the contraction 
matrix x, (^) and (|36D . 

The elements of the decomposition matrix are expressed as single finite sums of 
quotients of Pochhammer symbols involving the space-time dimension n. Whenever one 
of the indices j or j' reaches its minimal value (i.e. ji = 0, j2 = 0, js = or js = 1), we 
have just a single term (eqs. (|27|) , (0) and (|33D ). In most cases, however, a sum of several 
terms cannot be avoided because the numerator of 4>jij2j3;j[j'2j'cf contains non-factorizable 
quadratic (or higher degree) factors, e.g.: 

+ 3n + 6 



a,l,2;l,l,2 



72(n - l)n(n + l)(n + 2){n + 4)(n + 6) ' 



(45) 



Therefore, it seems unlikely that more simple general formulae for ^ju^isyiiaia ^^^^ 
and (B2D can be found. 



Numerous checks on the results were performed by computer algebra ||21|. We also 



verified, for a large number of cases, that the reduction formula in ref. |15| is equivalent 



to ours. The reduction algorithm in is of a different kind, because it is based on 
integration by parts and depends on the particular form of the propagators in the 
diagrams. 

Let us now briefiy discuss the application of these formulae. As an example, we 
consider the momentum expansion of a massive two-loop three-point function (see e.g. in 

|15|). Here, one can choose between two methods. The first is to expand all propagators 
in powers of the external momenta, yielding a collection of vacuum integrals with tensor 
numerators. Then, all the tensor integrals are reduced to scalar integrals times invariant 
tensor structures {ji, j2, is}- After that, the integrals can be calculated by expressing the 
numerators in terms of inverse propagators and cancelling them against the denominators. 
The final step is to contract the tensor structures with the external vectors (momenta or 
polarization vectors), which gives us a polynomial in scalar products of those vectors. If 
two or more of the external vectors are identical, each term is accompanied by a symmetry 
factor which, in general, can be written as a multiple sum of combinations of factorials, 
or computed by a recursive algorithm. 

The second method is to decompose the three-point function into scalar form factors, 
which are then expanded in the external invariants by using d'Alembertian operators. It 
takes some work to apply the d'Alembertians to the integrand, but in return, we no longer 
need to contract any tensor structures with external vectors. After setting the external 
momenta to zero, we get scalar vacuum integrals, which are evaluated in the same way 
as before. For given Ni and N2, we end up with three sums: one over powers of the 
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external momenta, one over powers of the d'Alembertians, and the internal sum in the 
representation of 0. As we have already mentioned, two of these sums can be absorbed 
in Gegenbauer polynomials, in which case only one extra sum remains. 
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